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ABSTRACT 

Context. AG Draconis produces a strong supersoft X-ray emission. The X-ray and optical/UV fluxes are in a strict anticorrelation 
throughout the active and quiescent phases. 

Aims. To identify the source of the X-ray emission and reveal the nature of the observed flux anticorrelation. 

Methods. The X-ray and UV observations with XMM-Newton, far-UV spectroscopy from FUSE, low- and high-resolution IUE 
spectra and optical/near-IR spectroscopic and/or photometric observations. Modeling the spectral energy distribution and broad wings 
of the Ovi /U032, /U038 and Heii/11640 lines by the electron-scattering during the maximum of the 2003 burst, the following 
transition and quiescent phase. 

Results. The X-ray-near-IR energy distribution at different levels of the star's brightness confirmed quantitatively the observed flux 
anticorrelation and showed that the optical bursts are associated to an increase of the nebular component of radiation. The profile- 
fitting analysis revealed a significant increase in the mean particle density around the hot star from ~ 2.6 x 10'° cnT 3 during quiescent 
phase to ~ 1.1 x 10 12 cirr 3 during the burst. 

Conclusions. The supersoft X-ray emission is produced by the white dwarf photosphere. The X-ray and far-UV fluxes make it 
possible to determine its temperature unambiguously. The supersoft X-ray-optical/UV flux anticorrelation is caused by the variable 
wind from the hot star. The enhanced hot star wind gives rise to the optical bursts by reprocessing high-energy photons from the 
Lyman continuum to the optical/UV. 
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1. Introduction 

Symbiotic stars are long-period (orbital periods are in order of 
years) interacting binaries consisting of a cool giant and a com- 
pact star, which is in most cases a white dwarf (WD), that ac- 
cretes from the giant's wind. This process generates a very hot 
(7h ~ 10 5 K) and luminous (Lh ~ 10 2 - 1O 4 L0) source of 
radiation, that ionizes a fraction of the neutral wind from the 
giant giving rise to nebular emission. As a result the spectrum 
of symbiotic stars consists of three basic components of radia- 
tion - two stellar and one nebular, (see e.g.lSeaauist et al Hll984t 
iKenvon & Webbink 1 1 1 984t iMtirset et ail I 1 S 



Skopal 2005). 



This situation represents the so-called quiescent phase, during 
which the symbiotic system releases its energy approximately at 
a constant rate and at a stable spectral energy distribution (SED). 
Sometimes, the symbiotic system changes its radiation signifi- 
cantly, brightens up by 1-3 magnitudes in the optical, shows sig- 
natures of a mass-outflow and c hanges its ionization structure for 
a few months to year s (see e.g. I Viotti et al. |[T984tlCorradi et al. I 
120031: I Skopal I I2005L and references therein). We name this 
stage as the active phase, and the corresponding brightening 
in the light curve is classified as the 'Z And-type' outburst. 
Occurrence of these outbursts is unpredi ctable and their nature i s 
so far poorly understood (e.g. Sect. 1 in Sokoloski et al. 2006). 
Investigation of this type of outbursts represents the key problem 
in the research of symbiotic stars. 

A G Dra is a yellow sym biotic binary comprising a K2III 
giant (Mti rset & Schmidll999t) and a WD accreting from the gi- 



ant's wind on a 549-day orbit (Fek el et al. 120 00). The light curve 
of AG Dra s hows numerous bursts with am plitude of 1-3 mag 
in U (Fig. D. lGonzalez-Riestra et al. I {1999) identified cool and 
hot outbursts differing in their Zanstra temperatures and the light 
curve profile. The former are more pronounced (AU ~ 3 mag) 
lasting for 1-2 years (e.g. 1981-83, 2006-07), while the latter 
are weaker (AU ~ 1 - 2 mag), single brightenings lasting for 
weeks to months (e.g. 1985, 1986, 20 03, see Fig. 1). Modeling 
the UV/IR continuum, Skopal (2005) found a significant con- 
tribution from the nebula in the near-UV/optical that strength- 
enes during outbursts. The symbiotic nebula in AG Dra is dense. 
Material supplied by the giant's wi nd into the b inary environ- 
ment at a rate of 3 x 10~ 7 Mq yr~ l dSkopal 112005b corresponds 
to number densities of ^ 1-10x1 8 ctrT 3 in between the bi- 
nary components and the neighbouring regions of similar dimen- 
sions. Therefore, to investigate properties of symbi otic nebula 
we can consider th e Case B for its radiation (see Oste rbrockl 
ll 7 9l9llKwokll2000l) . 

AG Dra is a halo binary system, placed at a galac tic lati- 
tude of +41° with a low reddening. Viott i et al. I d 19831) found 
£ B -v = 0.06 + 0.03 by fitting t he 2 2 00 A depression on the IUE 
spectra. Mikolajewska e t al. I d 19951) re-analyzed a large num- 
ber of IUE spectra and suggested £b-v = 0.05 + 0.01, while 
Birri el et al. I {2000) derived the extinction value of £b-v = 
0.08 + 0.01 mag by fitting the UV spectra obtained by the 
Hopkins Ultraviolet Telescope. These quantities are comparable 
to the total reddening through the Galaxy towards AG Dra (0.07- 
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Table 1. Log of the used spectroscopic observations 



Date 


Julian date 


Stage* 


Re 


'ion 


Observatory 


* exp 


YYYY/MM/DD 


JD24... 




[nm] 




M 


2003/10/06 


52919.5 


HB 


330- 


1020 


Loiano 


660 


2003/10/10 


52923.7 


HB 


212 


XMM-OM 


1200 


2003/10/10 


52923.7 


HB 


1.93 


-6.92 


XMM-Newton 


17121+ 


1985/03/13 


46137.8 


HB 


115 


-335 


IUE 


600 


1995/07/28 


49927.5 


HB 


Hen 


164.0 


IUE 


5400 


2003/11/14 


52958.5 


T 


100- 


108.2 


FUSE 


4863 


2003/11/19 


52963.7 


T 


1.93 


-6.92 


XMM-Newton 


9217 


2003/11/19 


52963.7 


T 


212 


- 291 


XMM-OM 


4400 


2004/06/15 


53172.0 


Q 


3.12 


-6.92 


XMM-Newton 


12354 


2004/06/15 


53172.0 


Q 


212 


- 291 


XMM-OM 


4400 


2004/06/24 


53181.0 


Q 


100- 


108.2 


FUSE 


10832 


1980/05/23 


44383.1 


Q 


115 


-335 


IUE 


2520 



* hot burst (HB), transition to quiescence (T), quiescent phase (Q), * no detection of supersoft X-rays 



0.08, iBurstein & Heiles 1 11982|) . which justifies that the mea- 
sured extinction to AG Dra is purely of the interst ellar nature. 
From t he Lya width on high-resolution IUE spectra. lViotti et alTI 
( 1983) estimated the H i column d ensity of log /Vh = 20.2 . Using 
the Einstein X-ray observations, lAderson et al. I (Il98ll) found 
Nh = 3 x 10 20 ctrr 2 . According to the relatio nship between A^h 
and £ B -v (AW£b-v=4.93 x 10 21 cm" 2 mag~ 1 . lDiplas & Savage I 
1994), both these parameters are consistent. 

The orbital inclination of the binary is rather low. There are 
no signs of eclipses either in the optica l, far-UV or X-ray re- 
gions (e.g. Gonzal ez-Riestra et al. 120081). Considering geom etry 
of the main sources of radiation. iMikolajewska et alTI (Il995[) es- 
timate d a system inclination i « 30° -45° and lSchmid & Schild I 
dl997l) . based on spectropolarimetric observations, suggested 
i = 60 + 8° .2. As a result we see the hot star rather from its pole 
than the orbital p lane, where a disk-like structur ed material can 
be expected (e.g. iMastrodemos & Morrisl ri998). Both the low 
interstellar absorption/extinction and the low orbital inclination 
suggest a high transmission of the interstellar medium (ISM) to 
soft X-rays, which is consistent with the fact that AG Dra is 
the brightest system at these energies among other symbiotics. 
Therefore, AG Dra has been frequently observed by X-ray satel- 
lites. 

Using the ROSAT observations iGreiner et aTl d 19971) first 
noted a remarkable decrease of the X-ray flux during the op- 
tical (1994 and 1995) maxima, while in the UV and the op- 
tical they indicated a large increase of the emission line and 
continuum fluxes. They ascribed this X-ray/UV flux anticorre- 
lation to a temperature decrease of the hot component. The X- 
ray history of AG Dra was re cently reviewed and discussed by 
IGonzalez-Riestra et al. I (120081) . They found that the anticorrela- 
tion between the X-ray and optical/UV emission appears to be 
a general feature of AG Dra radiation and is independent of the 
type and strength of the outburst. They suggested that during 
outbursts the WD radiation increases, but is strongly absorbed 
by the circumstellar ionized gas. 

As the effect of the flux anticorrelation is related to tran- 
sitions between active and quiescent phases, its understanding 
thus can aid us in revealing the nature of the Z And-type out- 
bursts. Accordingly, in this paper we investigate the origin of 
the observed anticorrelation between the supersoft X-ray and the 
optical/UV fluxes for the case of AG Dra. For this purpose we 
model its composite continuum within the X-rays - near-IR do- 
main at different levels of the activity. In Sect. 2 we summarize 
and describe nearly-simultaneous observations we used to model 



the continuum. In Sect. 3 we describe our analysis and present 
the results. Their discussion and summary with conclusions are 
found in Sects. 4 and 5, respectively. 



2. Observations 

2.1. Sources of the data and their timing 

For the aim of this paper we selected observations from the 
supersoft X-ray to the near-IR, taken at different levels of the 
AG Dra activity. We analyze the archival supersoft X-ray data 
made by the X-ray Multi-Mirror Telescope (XMM-Newton), in- 
cluding its Optical Monitor (XMM-OM) for the near-UV fluxes, 
the far-UV spectra (988 - 1082 A) made with the LiFIA channel 
of the Far Ultraviolet Spectroscopic Explorer (FUSE), the ultra- 
violet low- and high-resolution spectra taken by the International 
Ultraviolet Explorer (IUE), the optical low-resolution spectrum 
from the Loiano observatory and the flux-points determined by 
the broad band optical UBVR c Ic and near-IR JHKLM pho- 
tometry. The XMM -Newton and Loiano ob s erva tions were de- 
scribed in detail by IGonzalez-Riestra et aTl d2008l) . The FUSE 



spectra were processed according to Skopal e t al. ( 120061) . 



During quiescent phase we composed the observed SED 
by the XMM-Newton and FUSE observations, complemented 
with the IUE spectrum (SWP9084/LWR07831) taken at a sim- 
ilar orbital position. During transition phase we selected the 
near-simultaneous XMM-Newton and FUSE observations and 
the photometric UBV flux-points from the descending branch 
of the 2003 burst (Fig. 1). To model the SED at the maximum 
of hot bursts we used the ultraviolet XMM-OM fluxes and the 
optical spectrum, both from the maximum of the 2003 burst. 
However, to estimate the relevant far-UV fluxes we had to use a 
non- simultaneous observation. We used well-exposed IUE spec- 
tra (SWP25443/LWP05513) taken around the maximum of the 
1985 hot burst. To match the XMM-OM fluxes we scaled the 
IUE spectrum by a factor of 1.3. This was possible, because hot 
bursts ar e similar in their profiles and co lours (see Fig. 1 here and 
Fig. 6 of IGonzalez-Riestra et al. Ifl 999L for UV colours). Their 
SED shows a dominant contribution from the nebula to the near- 
UV/U region, while t he hot stellar source dominates the far-UV 
region (see Fig. 13 of lSkopalll2005l) . 

For the profile-fitting analysis we used the O vi /11032, /U038 
doublet on the FUSE spectra, available from the quiescent and 
transition phases. For the optical maximum we analyzed the 
He ii A 1640 line exposed on the high-resolution IUE spectra 
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Fig. 1. The U and V light cur ves of AG Dra from 1977. They are characterized with a series of outbursts with multiple maxima. 
The data were summarized by Skopa l et al. I d2007l) . The right panel shows a detail around the 2003 hot burst with timing of the 
XMM-Newton (long thin lines) and the FUSE (short thick bars) observations. During the optical maxima, the soft X-ray emission 
was not detectable (long dashed lines). 



(SWP55372 and SWP55373) at the maximum of the 1995 hot 
burst. The spectra were calibrated with the aid of their low- 
resolution counterparts. 

Relevant observations were dereddened with £b-v - 0.08 
and resulting parameters were sc aled to a distance of 1.1 kpc 
dBirriel et al. 12 000: Skopal 2005|). Their log and plots are given 
in Table 1 and Fig. 2, respectively. 

2.2. Derivation of the X-ray fluxes 

In this paper we analyze o bservations of AG Dra made by 
XMM-Newton dJansen et al. 11200 ll) . performed during the 2003 
hot burst (at its optical maximum and transition to quies- 
cence) and during quiescent phase (Fig. 1). Details of the ob- 
servations and the data reduc tion were already described by 
IGonzalez-Riestra et al. I d2008l) . To estimate the X-ray fluxes we 
fitted the spectra from the EPIC-pn instrument with the XSPEC 
software package. According to a very hi gh temperature of the 
hot stellar source in AG Dra (> 10 s K, [K enyon & Webbink 
[l984l:lMurset et al. ll99U[Greiner et al. ll997HSkopall2005l) we 
can assume that the supersoft X-ray fluxes are emitted by the 
WD photosphere. Therefore, we considered a simple tempera- 
ture blackbody model in fitting our 0.1 8-0.4 keV data. The re- 
sulting fit to the spectrum from quiescence corresponded to the 
blackbody radiation absorbed with the hydrogen column den- 
sity Nu ~ 3 x 10 20 cm~ 2 and temperature of ~14.5eV as in the 
case of the ROSAT PSPC observations from q uiescence (see 
iGreiner et all 1 19971: IGonzalez-Riestra et al. IfeOQSl) . During the 
transition, the shape of the spectrum was similar, but the lower 
signal to noise ratio prevented a meaningful fit. In this case we 
assumed the same parameters as for the quiescence spectrum, 
and thus fitted only the normalization factor. Fluxes were derived 
from these 'unfolded' models and the corresponding data/model 
ratio. This approach implies that the fluxes obtained by this way 
are, to a certain extent, model dependent. However, we are con- 
fident that the spectrum can be well represented by a simple ab- 
sorbed blackbody model, which is also supported by the data- 
minus-model residuals that are small, flat and within the errors 
in the range of 40-7 A. Here we refer the reader to the paper of 
iNowak et al. I (2002), who discuss other approaches to obtaining 
flux points from X-ray data and the effects connected with. To 



assess influence of the hydrogen column density and the temper- 
ature in the derived fluxes we fitted the data with a few hundred 
of models covering a wide range of these parameters. Taking the 
30 best fit cases (in terms of^ 2 ), we found that the fluxes derived 
from these models differ from the resulting ones by less than 
15% in the range of 30-65 A. Uncertainties of the X-ray fluxes 
derived in this way satisfy the purposes of this work. Further, 
we complemented the X-ray fluxes with those from the FUSE 
spectra to determine final parameters of the radiation produced 
by the WD photosphere (see below, Sect. 3.1.1.). 

3. Analysis and results 

3.1. Modeling the composite spectrum 

The continuum spectrum of symbiotic stars is composed of three 
basic components of radiation - two stellar, F^(A) and F „(A), 
from the hot star and the cool giant, respectively, and one neb- 
ular, F n (A), from the ionized gas in the system (Sect. 1). Their 
superposition then yields the observed flux as 

F(A) = F g (A) + F h (A) + F n (A). (1) 

To achieve our aim we need to determine physical parameters 
of the hot stellar and the nebular component of radiation during 
different stages of the activity. We will use the method of dis- 
entangling t he composite sp ectrum of S-type symbiotic stars as 
proposed by Skopal] d2005l) . Observations from the soft X-rays 
to the near-IR cover well all the energy domains, within which 
the individual components of radiation have a dominant contri- 
bution. This allows us to model them separately. In the following 
three sections we describe our approximations of these compo- 
nents of radiation. 

3.1 .1 . The hot star continuum 

According to [Greiner et all d 19971) and IViotti et al. I d2005h 
the soft X-ray radiation in AG Dra could be associated with 
the hot star photosphere. However, the soft X-rays from cos- 
mic sources are significantly attenuated by absorptions in the 
ISM, whi ch requires relevant correction before their interpreta - 
tion (e.g. lCruddace et al. 1 1974 1 Wilms. Allen & McCrav 12 000). 
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Assuming that no emission occurs on the path between the X-ray 
source and the observer, we can use the simplest solution of the 
radiative transfer equation to correct the observed X-ray flux, 
F° b *(A), for absorptions as 



F° bs (A) = F x (A)e- 



(2) 



where F X (A) is the X-ray spectrum as emitted by the source and 
t,i is the optical depth along the line of sight. It is given by 
the absorption cross-section per atom, <Ti(A), of the element, i, 
and the total number of atoms on the line of sight, Thus, 
ta = Ei nm(X) = E; ai<Ti(A)Nn = o- x (A)N H , where a t is the 
relative abundance, Nu is the total hydrogen column density 
[cm -2 ] and cr x (A) [cm 2 ] is the total cross-section of the mate- 
rial o n the line of sight per hydrogen atom in the X-ray domain 
(e.g. lCruddace et al. Ill 974 . 

In addition to the X-ray fluxes from XMM-Newton, we use 
the far-UV continuum fluxes between Al 188 and /U000 A, made 
nearly-simultaneously with the FUSE satellite (see Fig. 1). As 
we analyze both the X-ray and the far-UV data, we consider at- 
tenuation of the light by bound-free absorptions in the X-ray do- 
main and by the interstellar extinction in the far-UV spectrum. 
According to previous models (see Sect. 2.2.), we approximate 
the hot stellar continuum by a blackbody radiation at a temper- 
ature TV As a result, and in the sense of Eq. (2), we model the 
hot star continuum by fitting the observed X-ray/far-UV fluxes, 
F° bs (A), with a function 



for,! < 912 A 



K W = \ , ( 3 ) 

\6lnB A (T h ) 10-°- 4s *' £ b-v for A > 912 A 



where the scaling factor 9^ = R e ® Id is given by its effective 
radius and the distance d, and thus represents the angular 
radius of the hot stellar source. The observ ed far-UV fluxes were 
corrected using the extinction curve kx of lCardelli et al. | fl 989) 
and th e ratio of total to selective extinction R — 3.1 (e.g. Wegner 
120031) . In addition to the attenuation effects included in Eq. (3), 
the neutral hydrogen on the line of sight causes a depression of 
the continuum around hydro gen lines of the L yman series due 
to the Rayleigh scattering (Islik eret al. I Il989h . By analogy to 
the first term of Eq. (3) we can express the Rayleigh attenuated 
continuum as 



F° bs (A) = %nB A (T h )e-' r * 



(X)N H 



(4) 



where cr R (/l) is the Rayleigh scattering cross-section for atomic 
hydrogen. On our FUSE spectra the effect is seen around the Ly- 
B line (Fig. 2), and can be used to test quantity of Nu obtained 
from the X-ray domain. 

3.1.2. The nebular continuum 

We approximate the nebular component of radiation in the 
UV/optical continuum by processes of recombination and ther- 
mal bremsstrahlung in the hydrogen plasma, radiating under 
conditions of the Case B. Validity of this simplification is sup- 
ported by the following arguments: (i) A strong nebular contin- 
uum is characterized with a simple and low electron temperature 
of T e ~ 2 x 10 4 K (Table 3), which suggests that the nebula 
is powered mainly by photoionization, i.e. ionizations by colli- 
sions are not important, (ii) There are no recognizable signatures 
of the Hen continuum in the IUE spectra (e.g. a jump in emis- 
sion at ~2050 A and/or a pronounced Paschen series of the He n 



recombination lines). Also it is not possible to separate contri- 
butions from the He i continuum, because o f its very similar re- 
comb ination coefficients to that of Hi (e.g. iBrown & Mathews I 
1 1970b and a small abundance. That is why we consider nebular 
emission from hydrogen only, (iii) Symbiotic nebulae are rela- 
tively very dense (Sect. 1). This makes the mean-free path of 
any diffuse Lyman continuum photons to be too short to escape 
the nebula, i.e. the nebula is optically thick in the Lyman contin- 
uum. This means that the ionizations caused by stellar radiation- 
field photons are balanced by recombinations to excited levels 
of H i, while the ground state in the recombination process can 
be ignored. Therefore we consider the Case B for the nebular 
radiation, (iv) Because of the high density, we can neglect con- 
tributions due to the two-photon emission. 

Accord ing to these s implifications and with the aid of 
Eq. (11) in lSkopall d2005l) . the F R (A) term in Eq. (1) here can 
be expressed as 



EM 
And 2 



e*(H, T e ), 



(5) 



where e / j(H, T e ) [ erg cm 3 s _1 A -1 ] is the volume emission co- 
efficient for hydrogen, which depends on t he electron tempera- 
ture r e and is a function of wavelength (e.g. lBrown & Math ews 
Il970h . We note that the observed nebular emission re presents 
only the optically thin part of the symbiotic nebula (Skopal 
2001). This implies that the measured EM represents a lower 
limit of that originally created by ionizations, which puts a lower 
limit to the required flux of ionizing photons (see Eq. (6) below). 



3.1 .3. Radiation from the giant 

For the stellar radiation from the giant in AG Dra we adopt the 
model SED according to Skopal (2005). The model is based on 
the optical VRI and the near-IR JHKLM flux points matched 
by a synthetic spectrum calculated for the effective temperature 
of 4300 K. This spectrum then defines the first term in Eq. (1), 

9.51 x 10~ 9 erg ctrT 2 s~', 



Fg(A). Its bolometric flux, 2 crT* s 
corresponds to the giant's radius R g = 33 (d/l.l kpc) Rq and the 
luminosity L g = 360 (d/l.l kpc) 2 Lq. 



3.2. Thomson-scattering wings of O vi /11032, /U038 and 
He N/11640 lines 

The aim of this section is to model the extended wings of the O vi 
and He n lines measured at similar levels of the optical bright- 
ness, at which the X-ray observations were carried out. By this 
wa y we support the resu lts obtained by model SEDs. 

Schmi d et al. 1 (119991) suggested that the broad wings of the 
Ovi /U032 and /11038 resonance lines could be explained by 
scattering of the Ovi photons by free electrons. The effect of 
this process is weak and wavelength independent, because of 
a very small and constant value of the Thomson cross-section, 
<t x = 6.652 x 10 -25 cm 2 . However, the densest portions of the 
symbiotic nebula (log(n e ) ~ 8 — 12 cirT 3 ) extending to a few AU 
could be sufficiently optically thick for the electron-scattering 
process. Qualitatively, the effect should be more significant dur- 
ing active phases, because of a surplus o f electrons from the in- 
creased hot star wind (e.g.(Skopal 2006). Especially, the strong 
emission lines of highly ionized elements that are formed in the 
densest part of the hot stellar wind represent the best candidates 
for a well measurable effect of the Thomson scattering. 

To calculate the electron- scattering wing p rofiles we adopted 
the procedure suggested by Munch (1950$ that assumes the 
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Table 2. Parameters of the profile-fitting analysis (Sect. 3.2, 
Fig. 2) 



Date 


Stage 


Line 




T e 


2 

'tred 












[K] 


[cm- 3 ] 


1995/07/28 


HB 


Hen 


0.32 


28 800 


0.9 


1.1 x 10 12 


2003/11/14 


T 


Ovi 


0.083 


28 500 


1.2 


1.7 x 10" 


2004/06/24 


Q 


Ovi 


0.061 


20 300 


0.8 


2.6 x 10 10 



electron scattering happens in the layer outside the line for- 
mation region, and that the electrons are segregated from 
the other opacity sources, which implies no change in the 
eq uivalent width of th e line. Further we used the expression 
of ICastor et al. I (Il970b for the resulting scattered line profile 
and the electron-sca t tering redistribution function derived by 
lHummer~& Mihalas] (1 19671) . This simplified appr oach was used 
by many authors, recently bv lYoung et al. I ([2005). 

In our profile-fitting analysis, we approximated the observed 
P Cygni-type of the O vi line profiles by two Gaussians and fitted 
their broad wings with 2 free parameters - the electron tempera- 
ture, T e , and the electron-scattering optical depth, r e . Removing 
emission/absorption features and bumps in the profiles we fit- 
ted the broad wings in the range of about +10 A around the line 
cores using 644, 1004 and 370 flux-points in the wing profile 
from quiescence, transition and burst stage, respectively. The re- 
sulting fits for all three cases we analyzed have small reduced 
Xred ~ 1 (Table 2, Fig. 2), which implies that mean residu- 
als are comparable to the errors of the modeled flux-points. 
Thus the model fits well the broad wings of the profile, which 
confirms their origin as due to the electron-scattering. A de- 
tailed description of our approach will be presented elsewhere 
dSekeras& Skopal 112009b . 

A comparison of our model SEDs and profiles with observa- 
tions is plotted in Fig. 2 and corresponding parameters are given 
in Tables 2 and 3. 

3.3. The SED-fitting analysis 

3.3.1 . Parameters of the hot stellar source 

To derive parameters of the hot stellar source in AG Dra, the 
following points are relevant. 

(i) We dereddened the FUSE fluxes and calculated the 
Rayleigh scattering effect around Ly-/3 (Eq. (4)) to estimate Nn- 
The scattering cros s -sectio n, <x R (/t), was calculated according to 
iNussbaumer et al.1 dl989f) . However, numerous and strong ab- 
sorptions in the FUSE spectrum did not allow us to determine 
unambiguously the Nh parameter. Therefore, we could only 
compare the Rayleigh attenuated blackbody radiation to the ob- 
served continuum for a reasonable range of Nh quantities. We 
found that the observed depression in the continuum around Ly- 
j6 constrains Nh to ~ 2 - 4 x 10 20 cm -2 that is consistent with 
the values suggested by different methods (Sect. 1). 

(ii) We calculated a grid of synthetic models of the func- 
tion (3) for reasonable ranges of the fitting parameters, 0^, Tj, 
and Nh- We compared models to both the observed X-ray and 
dereddened UV fluxes, and used the x 2 statistics to evaluate 
the fit. To correct the observed X-ray fluxes for absorptions we 
used the tba bs absorption model (Tubingen- Boulder absorption 
ISM model. IWilms. Allen & McCrav I feOOO). During the quies- 
cent phase we fitted 1 1 X-ray fluxes from 31 to 69 A and three 
far-UV flux-points around 1000 A (bottom left panel of Fig. 2). 



We estimated errors in the X-ray fluxes to be of 2 — 17% of their 
mid values, while for the FUSE spectra we adopted errors as 
large as 5 — 10% of the measured continuum. The best-fit-model 
and the flux-point errors yielded^ 2 = 34.6 and^ 2 ed = 3.1 for 
11 degrees of freedom (dof). As the model fits well the data, 
the somewhat larger x 2 ed reflects rather small values of our error 
estimates. During the transition stage we fitted 7 X-ray fluxes 
between 43 and 69 A and three far-UV flux-points. We omitted 
fluxes around 20-30 A, because of their different nature. They 
could be originated by shocks in the nebula. Flux uncertainties 
we used, were estimated to be in the range of 7 — 22%. The re- 
sulting fit has a small x 2 ed = 0.6 for 7 dof, which suggests that 
the model fits the data very closely, well within their uncertain- 
ties. Corresponding best-fit-model parameters and their deriva- 
tives (the effective radius and the luminosity of the hot stellar 
source) are introduced in Table 3. 

(iii) During the bursts, when the supersoft X-ray emission is 
absorbed entirely, it is not possible to determine unambiguously 
the hot star temperature for Zj, ^ 10 5 K, because of too small 
range of the far-UV fluxes available. Nevertheless, we can es- 
timate a lower limit of the temperature of the ionizing source, 
r™ n , at which the hot star radiation, scaled to the far-UV fluxes, 
is just capable of producing the observed EM, i.e. the total num- 
ber of ionizing photons just balances the total number of re- 
combinations. According to simplifications for the plasma ra- 
diation as introduced in Sect. 3.1.2., we solve equilibrium equa- 
tion (Eq. 6 below) for the temperature and the scaling factor 
9^, which determine the flux of ionizing photons, with the input 
parameters of the measured nebula, EM and T e , The resulting 
temperature then corresponds to T™ n , at which the hot star ra- 
diatio n gives rise to the observed nebular emission (see[Skopal 
12001 in detail). 

3.3.2. Parameters of the nebula 

According to Eq. (1), the nebular component of radiation can 
be obtained by subtracting the stellar contributions from the ob- 
served spectrum, i.e. F n (A) = F(A)-Fh(A)-F g (A). First, we esti- 
mated representative continuum fluxes at about 20 wavelengths 
between 1250 and 3300 A and complemented them with pho- 
tometric U and XMM-OM fluxes, if available. Uncertainties 
of the well exposed spectra from the IUE archive are between 
5 and 10% of the measured values. Errors in the XMM-OM 
fluxes were estimated to only a few percents (see Table 2 in 
iGonzalez-Riestra et al. 1 120081) . Uncertainties of classical pho- 
tometric measurements can be assumed to be less than 10%. 
We corrected corresponding fluxes for t he influence o f emis- 
sion lines using our Loiano spectra (see ISkopal 1120071 in de- 
tail). Second, we subtracted the contribution of the WD and the 
giant and fitted the function F n (A) given by Eq. (5) to the cor- 
rected flux-points to determine the EM and T e of the symbiotic 
nebula. Resulting fits have x 2 ed ^ 1 that proves assumptions of 
our simplified model of the nebula (Sect 3.1.2.). Corresponding 
plots and parameters are in Fig. 2 and Table 3, respectively. 

4. Discussion 

4.1. Model SEDs and the flux anticorrelation 

Model SEDs show that the light variations in the optic a/near-UV, 
as measured at different levels of activity, are caused exclusively 
by the variable nebular continuum (Fig. 2, Table 3). Contribution 
from the hot stellar object can be neglected within this domain, 
and that from the giant can be assumed to be constant. 
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Table 3. Parameters of the SED-fitting analysis (see Sect. 3.3, Fig. 2). 
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3.1xl0 46 
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21500 


2.2xl0 59 
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Q 


2.82xl0 20 


164400 


6.3xl0~ 13 


0.031 


630 


3.4xl0 46 


3.1 


21000 


1.3xl0 59 
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The nebular component of radiation represents a frac- 
tion of the hot stellar radiation transformed by the ioniza- 
tion/recombination events throughout the symbiotic nebula. The 
result of this process depends on the number of ionizing photons 
(L p h [s -1 ]) produced by the hot star and the number of parti- 
cles on their path that are subject to ionization. In the case of 
the hydrogen plasma (see Sect. 3.1.2.) characterized with one r e 
and a mean particle concentration, «, the equilibrium condition 
between Lph photons and the number of recombinations can be 
expressed as 



Lp h (H) = f n e n + a B (H, T e ) dV = a B (H, T e ) h 2 V 
Jv 



(6) 



= a B (H,L e )LM, 



where as(H, T e ) [cm^s -1 ] is the recombination coefficient to 
all but the ground state of hydrogen (i.e. the Case B). For pa- 
rameters of AG Dra during quiescence, the equilibrium condi- 
tion is satisfied at the locus of points in directions from the hot 
to/around the cool star, where the flux of L p h photons are bal- 
anced by the flux of the neutral atoms of hydrogen in the wind 
from the giant. Throughout the remainder part of the nebula 
L p h > a^EM, which means that a fraction of L p h photons es- 
capes the nebula without being tran sformed to the nebular radi- 
ation (see Appendix B and Fig. 3 infSkopal 2001). Under these 
conditions, new particles injected into such the nebula will con- 
sume the surplus of ionizing photons, what consequently will in- 
crease the nebular emission observed during the bursts. This sit- 
uation can naturally be explained by a strengthened wind from 
the WD, that also induces an increase of the optical depth in 
the continuum in the direction of the WD, what we indicate by 
its larger effective radius in comparison with values from quies- 
cence (Table 3). 

According to Eq. (2) and the relatively large values of cr x in 
the supersoft X-ray domain, a small increase of A^h in between 
the source and the observer produces a significant attenuation of 
the F X (A) fluxes. For example, comparison of our XMM-Newton 
observations from the transition and quiescent phases shows that 
the increase of the A^h value by a factor of only 1.16 (Table 3) 
produces attenuation of the observed fluxes at ~ 0.18keV by a 
factor of ~ 5 (see Fig. 2). At the maximum of the 2003 burst, 
the 0.19 - 0.4 keV emission was not detectable, in spite of the 
increase of both the luminosity and the temperature of the hot 
source (Table 3, Fig. 2). The high limiting quantities of Lh and 
7h are constrained by the minimum flux of photons capable of 
ionizing hydrogen, L p h(H) = 8.9 x 10 46 s , that is required to 
give rise the observed high amount of the emission measure, 
EM = 6.2 x 10 5 V/l.lkpc) 2 cirT 3 (see Eq. (6)). We found that 
the value of A^h ~ 2.5 x 10 21 cirT 2 is sufficient to damp down the 
model fluxes (w S.SxlO^ergcm^s- 1 within the 0.1 9-0.4 keV 
range) below a detection level of around 10~ 16 erg cirT 2 s _I A -1 
for the exposure time used by XMM-Newton on 2003/10/10. 

According to the model SEDs we ascribe the observed flux 
anticorrelation to variable wind from the hot star. The wind 



particles enrich the plasma surrounding the WD's photosphere, 
which increases the number of both the bound-free absorptions 
(parametrized by A^h in the model) and the free-bound emis- 
sions. The former attenuates the supersoft X-ray fluxes, while 
the latter increases the nebular emission. 



4.2. Wing profiles and the electron density around the WD 

Enhancement of the particle concentration around the hot star is 
supported independently by our modeling the electron-scattering 
wing profile of the O vi and He n lines. The model parameter, the 
optical depth r e of the electron-scattering layer with the thick- 
ness r is related to its mean electron concentration n e as 



<Tjn e r. 



(7) 



In spite of the model simplification (Sect. 3.2) we approximate 
the thickness r of the scattering layer with the radius of the line 
emitting zone. Having its radius then allow us to estimate the h e 
quantity with the aid of Eq. (7). For example, the /?ovi radius 
can be obtained from the equation of the equilibrium between 
the flux of Lph(Ov) photons, capable of ionizing Ov ions, and 
the number of recombinations by O vi ions. By analogy with 
Eq. (6), and assuming spherically symmetric nebular medium 
around the central ionizing source, the equilibrium condition can 
be expressed as 



4-7T 

L ph (Ov) = —Rl, 



A(Ovi)0.83« 2 a B (Ovi), 



(8) 



where the factor 0.83 is the ratio of protons to electrons, A(0 vi) 
is the abundance of Ovi ions, aB(Ovi) is the recombination 
coefficient of a free electron with the O vi ion for the Case B. 
Substituting the radius of the O vi zone from Eq. (7) to (8), the 
electron concentration can be expressed as 



4-7T / T e 

Re = 0.83— — 
3 \cr T 



A(Ovi) 



«b(Ovi) 
WOv)' 



(9) 



The Lph(Ov) quantity as a functi on of the ionizing s ource tem- 
perature is showed in Fig. B.l of lSkopal et al. I (|2006). 

Our model parameters, Lh and 7h from quiescent and 
transition phase (Table 3) correspond to L p h(Ov) = 4.4 and 
1.7x10 s -1 , respectively. Assuming that all oxygen atoms are 
ionized to O vi within the zone, i.e. A(Ovi) = A(O) = 4.6 x 



10~ 4 dAsplundetal.l l2004. a B (Ovi) = 9.5 x lO ^cm^s 1 



dGurzadvan 1 1 1 997l) ~and T e from Table 2, yield n e = 2.6 X 
10 10 cm" 3 (R vi=50tf ), and 1 .7x 10 1 1 cirT 3 (R 0vl ±l 1 Rq) dur- 
ing quiescence (2004/06/15-24) and the transition from the burst 
(2003/1 1/14-19), respectively. 

During the optical maximum of the 2003 hot burst no FUSE 
observation was available. Instead, we modeled the Heii/11640 
line from the maximum of the 1995 hot burst (Table 1), because 
this line is in major part also created at a vicinity of the hot star. 
In contrast to the quiescent phase, an extended wing profile satis- 
fying the electron-scattering broadening developed at the bottom 
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of its intense emission core (~ 2.1 x 10 10 erg cm 2 s 1 A ') at 
the optical maximum (see Fig. 2, top right). A relatively high 
optical depth, r e = 0.32 (Table 2) L ph (Hen) = 2.7 x 10 46 s~' , 
or B (Heil) = 7.5 x lO^cm+V dNussbaumer & Vogel |[l987h 
and A(Hen) = 0.1, yield n e = 1.1 x 10 12 cnT 3 (R Heu =6.QR Q ). 
The presence o f the broad w i ngs of the He n A 1640 line was first 
pointed out by IViotti et aTl (Il983l) on the IUE high-resolution 
spectra from the 1981 outburst of AG Dra. Simultaneously, they 
found signatures of the hot star wind suggested by the P Cygni 
profiles of the N v /11238, /11242 resonance lines. 

4.3. Flux anticorrelation and the nature of bursts 

Results of the previous two subsections are mutually comple- 
mentary. Both the model SEDs and the profile-fitting analysis 
indicate a significant increase of the particle density around the 
WD during the bursts of AG Dra. According to other i ndepen- 
dent analyzes (see e.g. the Ho- method employed by Skopal 
2006), such the increase in the particle density is due to enhanced 
wind from the hot star into the surrounding particle-bounded 
nebula. This also explains the nature of this type of the opti- 
cal bursts, because new emitters will convert the excess of hy- 
drogen ionizing photons in the nebula into the nebular radiation 
that dominates the optical/near-UV. By return, this mechanism 
explains the origin of the X-ray-optical/UV flux anticorrelation 
as the result of variations in the wind from the hot star during 
different levels of the star's activity. 

Having the result of the reprocessing mechanism, i.e. the 
parameters of the nebular emission and/or of the wing profiles 
during bursts and quiescence, we can determine the mass loss 
rate from the accretor. Generally, the increase of the mass-loss 
rate via the wind can result from an increase in the mass accre- 
tion rate due to an accretion-disk instability. In such the case 
the enhanced wind provides an important mechanism for re- 
moval of angular m omentum of rapidly accreted material ont o 
the WD surface (e.g. lDuschl |[l9lllWarnerl[T9^invio1[T997l) . 
From this point of view, the origin of the inverse flux correlation 
suggests directions for further investigation of the nature of the 
Z And-type outbursts. 

5. Summary and conclusions 

In this paper we investigated the origin of the supersoft X-ray — 
optical/UV anticorrelation observed during different levels of ac- 
tivity of the symbiotic binary AG Dra. 

We modeled the composite continuum from the supersoft X- 
rays to near-IR during the maximum of the 2003 burst, the fol- 
lowing transition to quiescence and a quiescent phase. We de- 
termined physical parameters of individual components of radi- 
ation (Fig. 2, Table 3). To support the results obtained by model 
SEDs we fitted the broad wings of the O vi ,11032, /H038 and the 
Hen /11640 lines by the Thomson-scattering process. Resulting 
profile fits and corresponding parameters are shown in Fig. 2 
and Table 2, respectively. Main results can be summarized as 
follows. 

1. During the quiescent phase the radiation of the hot stellar 
source can be reproduced by that of a black body with a 
radius Rh = O.O31(c//l.lkpc)7?0 and radiating at the tem- 
perature 7h = 164 400 K, which yields the luminosity Lj, = 
630(a7 1. lkpc) 2 Lq. The X-ray emission was attenuated by 
absorptions corresponding to the neutral hydrogen column 
density N H = 2.82 x 10 20 cirT 2 for the ISM abundances 
that is equivalent to its interstellar value (cf. Sect. 1). The 



emission measure of the nebular component of radiation was 
EM = 1.3 x 10 59 (ii/l.lkpc) 2 cm" 3 . Fitting parameters for 
the electron-scattering wings of the O vi doublet (Table 2) 
and the number of L p h(Ov) photons from the model SED 
(Sect. 4.2) correspond to the mean electron concentration 
around the hot star, n e ~ 2.6 x 10 10 ctrr 3 . 

2. During the transition phase these parameters changed to 
R h = 0.036(d/l.lkpc)/? o , T h = 146 200 K, L h = 
53O(ii/l.lkpc) 2 L0, and the absorption of the X-ray fluxes 
corresponded to Nu = 3.27 x 10 20 cirT 2 . The larger quantity 
of the Nu parameter reflects an increase in the bound-free 
absorptions due to an enhanced wind from the hot star. As 
a result the X-ray fluxes decreased relatively to their values 
from quiescence. In contrast, the nebular emission increased 
to EM = 2.2 x 10 5 V/l.lkpc) 2 cm- 3 . Simultaneously, flux 
of the broad Ovi wings was by a factor of ~ 2 larger 
than during quiescence and the model parameters (Tables 2 
and 3, Eq. (9)) yield a significantly larger density of n e ~ 
1.7 x 10 11 crrT 3 around the hot star, within the O vi zone. 

3. During the burst, the high limiting quantities, 7^ = 
180000K and L h = 1760(^/1. lkpc) 2 Lq, are required 
to produce the observed large amount of EM = 6.2 x 
10 59 (c//l.lkpc) 2 cm -3 . The negative detection of the super- 
soft X-ray emission constraints a significant absorption ef- 
fect that can be parametrized with Nu ~ 2.5 x 10 21 cirT 2 . The 
strong and extended wings of the Hen/11640 line (Fig. 2) 
imply a large value of n e ~ 1.1 x 10 12 cirT 3 (Eq. (9)). 

Based on these results we formulate the following conclusions: 

(i) The model SED demonstrated that the supersoft X-ray 
emission is produced by the WD photosphere. The X-ray and 
far-UV fluxes allow us to determine its temperature unambigu- 
ously. 

(ii) The WD's continuum spectrum is modified by the cir- 
cumstellar and interstellar material in the line of sight. 

(iii) We found that the source of the opacity, causing the ob- 
served anticorrelation between the X-ray and optical/UV fluxes, 
can be associated with the hot star wind , which enhan ces during 
active phases of symbiotic binaries (see Skopal 2006). 

(iv) The higher mass loss rate increases the particle density 
at the vicinity of the WD. This event increases the number of 
bound-free absorptions in the line of sight, which leads to a sig- 
nificant attenuation of the supersoft X-ray photons and, conse- 
quently, the free-bound transitions under the Case B increases 
the nebular emission that dominates the optical/near-UV. 

(v) The origin of the X-ray-optical/UV flux anticorrela- 
tion explains the nature of bursts in AG Dra by reprocess- 
ing high-energy photons into the optical through the ioniza- 
tion/recombination events. Understanding the inverse relation- 
ship between optical and X-ray fluxes represents an important 
ingredient in the investigation of the Z And-type outbursts. 

Acknowledgements. This work is in part based on observations obtained with 
XMM-Newton, an ESA science mission with instruments and contributions di- 
rectly funded by ESA Member States and NASA. The far ultraviolet data pre- 
sented in this paper were obtained from the Multimission Archive at the Space 
Telescope Science Institute (MAST). They were made with the NASA-CNES- 
CSA Far Ultraviolet Spectroscopic Explorer. FUSE was operated for NASA by 
the Johns Hopkins University under NASA contract NAS5-32985. The optical 
spectra were obtained from the archive of the Loiano Station of the Bologna 
Atsronomical Observatory. This research was in part supported by a grant of 
the Slovak Academy of Sciences No. 2/7010/27. The authors are grateful to the 
anonymous referee for critical comments and constructive suggestions. 

References 

Anderson, C. M., Cassinelli, J. P., & Sanders, W. T. 1981, ApJ, 247, L127 



8 



A. Skopal et al.: The origin of the supersoft X-ray-optical/UV flux anticorrelation in the symbiotic binary AG Draconis 



Asplund, M., Grevesse, N., Sauval, A. J., et al. 2004, A&A, 417, 751 
Brown, R. L., & Mathews, W. G. 1970, ApJ, 160, 939 
Burstein, D., & Heiles, C. 1982, AJ, 87, 1165 

Birriel, J.J., Espey, B.R., & Schutle-Ladbeck, R.E. 2000, ApJ, 545, 1020 

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 

Castor, J. I., Smith, L. E, & van Blerkom, D. 1970, ApJ, 159, 1119 

Corradi, R. L. M., Mikolajewska, J., & Mahoney, T. J. 2003, Symbiotic Stars 

Probing Stellar Evolution, ASP Conf. Ser. 303 (San Francisco: ASP) 
Cruddace, R., Paresce, E, Bowyer, S., & Lampton, M. 1974, ApJ, 187, 497 
Diplas, A., & Savage, B. D. 1994, ApJ, 427, 274 
Duschl, W. J., A&A, 163, 56 

Fekel, F. C, Hinkle, K. H., Joyce, R. R., & Skrutskie, M. 2000, AJ, 120, 3255 
Gonzalez-Riestra, R., Viotti, R., Iijima, T., & Greiner, J. 1999, A&A, 347, 478 
Gonzalez-Riestra, R., Viotti, R. F., Iijima, T., Rossi, C, Montagni, E, Bernabei, 

S., Frasca, A., Skopal, A. 2008, A&A, 481, 725 
Greiner, J. Bickert, K., Luthardt, R., Viotti, R., Altamore, A. Gonzalez-Riestra, 

R., Stencel, R. E. 1997, A&A, 322, 576 
Gurzadyan, GA. 1997, The Physics and Dynamics of Planetary Nebulae (Berlin: 

Springer- Verlag), p. 50 
Hummer, D. G., & Mihalas, D. 1967, ApJ, 150, L57 
Isliker, H., Nussbaumer, H., & Vogel, M. 1989, A&A, 219, 271 
Jansen, E, Lumb, D., Altieri, B., et al. 2001, A&A, 365, LI 
Kenyon, S. J., & Webbink, R. F. 1984, ApJ, 279, 252 

Kwok, S. 2000, The Origin and Evolution of Planetary Nebulae (Cambridge: 

Cambridge University Press), p. 13 
Livio, M. 1997, in ASP Conf. Ser. 121, Accretion Phenomena and Related 

Outflows, ed. D. T. Wickramashinge , G. V. Bicknell, & L. Ferrario, (San 

Francisco: ASP), 845 
Mastrodemos, N., & Morris, M. 1998, ApJ, 497, 303 

Mikolajewska, J., Kenyon, S. J., Mikolajewski, M, Garcia, M. R., & Polidan, R. 

S. 1995, AJ, 109, 1289 
Munch, G. 1950, ApJ, 112, 266 

Miirset, U., Nussbaumer, H., Schmid, H. M., & Vogel, M. 1991, A&A, 248, 458 

Miirset, U., & Schmid, H. M. 1999, A&AS, 137, 473 

Nowak, M. A., Wilms, J., & Dove, J. B. 2002, MNRAS, 332, 856 

Nussbaumer, H., & Vogel, M. 1987, A&A, 182, 51 

Nussbaumer, H., Schmid, H. M., & Vogel, M. 1989, A&A, 211, L27 

Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic 

Nuclei (San Francisco: W. H. Freeman and Company Press) 
Schmid, H. M., & Schild, H. 1997, A&A, 321, 791 
Schmid, H. M., Krautter, J., Appenzeller, I., et al. 1999, A&A, 348, 950 
Seaquist, E. R., Taylor, A. R., & Button, S. 1984, ApJ, 284, 202 
Sekeras, M., & Skopal, A. 2009, (in preparation) 
Skopal, A. 2001, A&A, 366, 157 
Skopal, A. 2005, A&A, 440, 995 
Skopal, A. 2006, A&A, 457, 1003 
Skopal, A. 2007, New. Astron., 12, 597 

Skopal, A., Vittone, A. A., Errico, L., et al. 2006, A&A, 453, 279 
Skopal, A., Vaflko, M., Pribulla, T. et al. 2007, Astron. Nachrichten, 328, 909 
Sokoloski, J. L., Kenyon, S. J., Espey, B. R., et al. 2006, ApJ, 636, 1002 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 

Viotti, R., Ricciardi, O., Ponz, D., Giangrande, A., Friedjung, M., Cassatella, A., 

Baratta, G. B., & Altamore, A. 1983, A&A, 1 19, 285 
Viotti, R., Altamore, A., Baratta, G. B., Cassatella, A., Friedjung, M. 1984, ApJ, 

283, 226 

Viotti, R. E, Gonzalez-Riestra, R., Iijima, T., et al. 2005, ApSS, 296, 435 
Warner, B. 1995, Cataclysmic variable stars (Cambridge, New York: Cambridge 

University Press) 
Wegner, W. 2003, AN, 324, 219 

Young, P. R., Dupree, A. K, Espey, B. R., Kenyon, S. J., & Ake, T. B. 2005, 
ApJ, 618, 891 



A. Skopal et al.: The origin of the supersoft X-ray-optical/UV flux anticorrelation in the symbiotic binary AG Draconis 



9 




l 1 1 1 1 1 1 1 r 




1.5 2 2.5 3 3.5 4 4.5 1020 1030 1040 

log(Wavelength) [A] Wavelength [A] 



Fig. 2. Left panels show a comparison of the observed and modeled SEDs of AG Dra during the hot burst (top), transition phase 
(middle) and quiescent phase (bottom) with corresponding residuals. Open/filled squares are the observed/corrected X-ray fluxes. 
They are in units of 1CT 12 erg cirT 2 s A -1 . Typical uncertainties of the flux-points are around of 10% (Sect. 3.3). Right panels 
compare the observed and modeled broad wings of the O vi /U032, /11038 doublet and the He n /U640 line at these stages of activity. 
Sample points with uncertainties are plotted in red. Timing of individual observations are given in Table 1. Models are described in 
Sect. 3 and their parameters are given in Tables 2 and 3. 



